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Summary 

Fluctuating pressure loads were measured along 
the top interior wall of a generic scramjet engine 
model during a series of performance tests conducted 
in the Combustion-Heated Scramjet Test Facility at 
the Langley Research Center (formerly the Langley 
Mach 4 Scramjet Test Facility). This facility can 
simulate Mach number 4 flight conditions at dynamic 
pressures that correspond to altitudes from 57000 
to 86000 ft. The 75-in-long model was installed in 
the test cabin of the facility with the inlet positioned 
at the tunnel nozzle exit plane. The nominal inlet 
Mach number was 3.5 at a total temperature and 
pressure of 1640°R and 92 psia, respectively. 

To minimize gage exposure to damaging heat 
flux from the high-enthalpy flow, 0.1 59- in-diameter 
piezoresist ive strain gages were recess mounted in 
existing calorimeter ports at nine locations along 
the interior of the model’s top upper wall. Recess 
distance ranged from 4 to 12 gage diameters. This 
technique provided a means to obtain meaningful 
data in the hostile flow environment. 

Fluctuating wall pressures were successfully mea- 
sured to 5 kHz at five locations upstream of the com- 
bustor region and at one location in the combustor 
section. However, not all gages were used for every 
run. The root mean square (rms) pressures ranged 
from 0.06 to 7.8 percent of the local dynamic pres- 
sure q e at the boundary layer edge. The rms pres- 
sures, mainly in the isolator section, greatly exceeded 
those reported in the literature for attached, turbu- 
lent boundary layers on flat plates, which generally 
tend to a limit of 1 percent q e . 

The overall fluctuating pressure levels in decibels 
(re 20 /xPa), over a frequency range 0 to 5 kHz, 
varied from 131 dB in the inlet region to 180 dB in 
the isolator-combustor region. These generally high 
levels were measured both before and after an inlet 
unstart. The main effect of the unstart was to extend 
the high-pressure levels further upstream from the 
combustor into the isolator section. 

A key finding was that combustion noise propa- 
gates upstream into the isolator section through the 
subsonic portion of the boundary layer. Combustion 
noise contributes to the total noise measured in the 
isolator section upstream of the combustor, where 
noise levels increased by 15 dB when the equivalence 
ratio (fuel-air ratio normalized by the stoichiometric 
ratio) was increased from 0.37 to 0.57. 

Transient pressure disturbances associated with 
thermally induced inlet unstarts were also measured. 
The unstart shock was found to propagate upstream 


at approximately 7 ft/sec, which suggests that tran- 
sient loading effects on the engine structure have 
little effect on structural damage compared with 
boundary layer combustion noise. 

Symbols 

A, B temperature-dependent parameters 

in gage transfer function, also 
empirical constants describing 
functional dependence of Strouhal 
number 

Ap y Bt derivatives of A and B with respect 

to temperature 

C empirical constant related to 

effective port area, 0.9 

c sound speed 

c e sound speed at boundary layer edge 

d cylindrical recess diameter 

E gage output voltage, instantaneous 

value 

E m mean value of gage output voltage 

/ frequency 

G power spectral density 

K v empirical constant (0.57) used in 

appendix C 

k poly tropic constant, also statistical 

degrees of freedom 

Z cylindrical recess length from port 

to gage diaphragm 

lij + 1 length of element (z, i 4* 1) 

M local Mach number 

M e local Mach number at boundary 

layer edge 

N number of sample points used in 

PSD analysis 

N ft number of data analysis blocks 

Nji e unit Reynolds number at boundary 

layer edge 

n cavity cross mode number, also nth 

sampling point 

p pressure 

P,Pm mean value of fluctuating pressure, 

also labelled p\ 

p fluctuating pressure 



■pin) 

pressure at nth sampling point 

Pit) 

pressure time history 

Pe 

mean static pressure at boundary 
layer edge 

Pg 

fluctuating pressure at gage face 

Pi 

fluctuating pressure at station i of 
transmission line 

Pirns 

rms pressure 

Ps 

mean static pressure at port 
entrance 

Ps 

fluctuating static pressure at port- 
entrance 

Pt 

mean total pressure 

Qe 

dynamic pressure at boundary layer 
edge 

R 

gas constant, 1716 ft 2 /sec 2 -°R 

R f 

recovery factor 

S p (u>) 

averaged power spectral estimates 

S p {lo) 

power spectral estimate of S p (lj) 

Sp(v) 

power spectral estimate for jth data 
block 

s 

cross-sectional area of port 

T 

fluctuating temperature 

Taw 

adiabatic wall temperature 

T e 

static temperature at boundary 
layer edge 

T t 

total temperature of gas 

T w 

wall temperature 

t 

time 

tB 

time length of data block 

^total 

time length of Nq contiguous data 
blocks 

U e 

gas velocity at boundary layer edge 

u 9 

fluctuating volume velocity at gage 
face 

Ur 

fluctuating volume velocity at 
station i of transmission line 

V 

cavity volume over gage diaphragm 

W s 

combined data window correction 
and calibration constant 


Z c tube characteristic impedance 

a empirical constant used in appen- 

dix C, 0.24 

<23 skewness 

Q4 kurtosis 

7 specific heat ratio 

6* boundary layer displacement 

thickness 

A increment 

fi\ mean value 

//3 third moment of pressure probabil- 

ity density about the mean 

/Z4 fourth moment of pressure probabil- 

ity density about the mean 

p cav gas density in recess cavity 

p e gas density at boundary layer edge 

p s gas density at wall 

a standard deviation, also gage 

diaphragm deflection factor 

a 2 variance 

cf) equivalence ratio, also propagation 

constant 

l o angular frequency 

oo free stream 

Abbreviations: 

ESP electronic scanning pressure system 

GHSE generic high-speed engine 

OASPL overall sound pressure level com- 

puted from variance 

(OASPL) psp overall sound pressure level com- 
puted from PSD 

probability density 

probability density function 

power spectral density 

absolute pressure in pounds per 
square inch 

root mean square 

run 6, batch 7 

run 15, batch 6 

run 50, batch 10 


PD 

PDF 

PSD 

PSIA 

rms 

R6B7 

R15B6 

R50B10 
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TDR transient data recorder 

TN tunnel nozzle 

Introduction 

The use of air-breathing propulsion systems for 
aerospace- plane vehicles requires a flight path opti- 
mized for maximum tolerable dynamic pressure to 
achieve high propulsive efficiency. Unfortunately, 
high dynamic pressures imply severe unsteady pres- 
sure loads on surface areas exposed to turbulent 
boundary layers. In localized areas of shock impinge- 
ment or flow separation associated with high sur- 
face curvature, unsteady pressure loads may be ex- 
acerbated such that overall levels may reach at least 
170 dB (re 20 /xPa), a level that may severely damage 
structures. As yet, no analytical models provide the 
structural design engineer with needed predictive ac- 
curacies for these unsteady loads. In particular, how 
are the characteristics of unsteady boundary layer 
loads on the interior of scramjet engines affected by 
interacting shock wave and boundary layers, sepa- 
rated flow over surfaces of high curvature, and the 
combustion process? 

Scramjet propulsion engineers seek to relate en- 
gine performance to unsteady loads measured at var- 
ious locations on the interior walls. A particular in- 
terest is the highly undesirable phenomenon known 
as “engine unstart, ” which is preceded by a destabi- 
lization of the inlet shock system. In the unstarted 
state, a bow shock forms in front of the engine inlet; 
the resulting flow pattern allows a large amount of 
spillage flow around the engine, with a consequent 
drag increase and combustor flame out. When the 
engine is operating at the highest possible level of 
heat release rate (i.e., near thermal choke), transient 
pressure disturbances from the combustion unsteadi- 
ness may initiate an unstart. High engine efficiency 
depends, in part, on the highest possible heat release 
rate consistent with stable operation. Thus, it is im- 
portant to quantify the unsteady loads in the engine 
interior before an engine unstart. 

Thus, two critical problems exist for designers of 
supersonic combustion ramjets (scramjets): control 
of the engine operating transients to prevent aerody- 
namic instabilities from initiating an engine unstart, 
and design of structures to withstand intense fluctu- 
ating pressure loads generated on the engine interior 
by aerodynamic and combustion processes. In his 
review of aerothermal problems associated with hy- 
personic flight, Holden (ref. 1) states that even the 
prediction of laminar-to-turbulent transition remains 
an important unresolved problem in hypersonic flow. 
Further, the complex flow fields around interacting 


shock waves and turbulent boundary layers make de- 
tailed analytical predictions of such regions almost 
impossible within the framework of the boundary 
layer equations. Therefore, an ability to predict fluc- 
tuating loads associated with these complicated flows 
must likely rely on a systematic measurement pro- 
gram to supply the necessary data base both for 
engine aerodynamics and component design. 

This work has a twofold purpose. The pri- 
mary purpose is to enhance and extend the exist- 
ing technology for fluctuating loads measurements 
over a broad frequency range underneath the bound- 
ary layers of high-speed and high-enthalpy flows. 
When exposed to such flows, the piezoresistive trans- 
ducer element typically used in conventional flush- 
mounted pressure gages sustains excessive thermal 
loads that cause gage sensitivity changes and even 
failure. Therefore, other measurement configurations 
must be found. Although active cooling has been a 
solution in some cases, this technique is too mechan- 
ically intrusive and inconvenient for most applica- 
tions. This paper documents an attempt to measure 
unsteady loads during performance tests of a generic 
scramjet engine model with high-frequency pressure 
gages recessed into the interior wall to reduce ther- 
mal flux to the gage diaphragm. The gages were 
installed in existing calorimeter ports. The model 
was designed by the Johns Hopkins Applied Physics 
Laboratories and tested in the Combustion-Heated 
Scramjet Test Facility at the Langley Research 
Center. 

The second purpose of the work is to present 
results of unsteady loads measurements in a scram- 
jet model. These data are relevant to aerody- 
namic performance and structural design. The re- 
sults include pressure time histories obtained during 
a thermally induced engine unstart, pressure spec- 
tra during steady-state operation, and the effects of 
increased heat release rate on fluctuating pressures 
in the boundary layer upstream of the combustor 
(isolator section). 

Experiment Configuration 

Test Facility 

The tests were conducted in the Combustion- 
Heated Scramjet Test Facility at the Langley Re- 
search Center (formerly the Langley Mach 4 Scramjet 
Test Facility). A schematic diagram of the facility 
is shown in figure 1. The facility is equipped with 
an air ejector to enable simulation of dynamic pres- 
sures from 500 to 1900 psfa (altitude from 86000 
to 57000 ft) at a Mach number of 4. Flight enthalpy 
is duplicated by the combustion of hydrogen in air 
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with oxygen replenishment to yield vitiated air that 
contains an oxygen volumetric content of 21 percent. 
The test gas from the burner is expanded through a 
converging-diverging nozzle to a Mach number of 3.5. 
This supersonic flow exhausts as a free jet from 
the 13-in. square nozzle exit into a cabin that houses 
the model. The gas passes through and around the 
model into an exhaust duct system connected to an 
annular air ejector that exhausts to the atmosphere. 
As shown in figure 1 , the facility is within a cell that 
has an ambient air intake tower at the upstream end 
and an exhaust tower at the downstream end. A 
more complete description of this facility, along with 
calibration results, is given in reference 2. 

Test Model 

A scramjet engine model, designed by the Johns 
Hopkins Applied Physics Laboratories, was installed 
in the facility cabin. Figure 2 shows the uninstalled 
model. The copper heat sink model (1200 lbm) 
was equipped with various spill doors and boundary 
layer bleed plates to condition the flow and to aid 
inlet starting. The copper mass provided, a heat 
sink for the intense thermal flux generated by the 
aerodynamic and combustion processes. In figure 3 
the model is shown installed in the facility with the 
spill window open. The flow is from right to left. Test 
results relating to the performance of this model as 
a propulsive device are described in reference 3 and 
are not discussed here. 

Important features of the model relevant to the 
present discussion are shown in the conceptual view 
of figure 4. The model plan view and side view 
(not to scale) are shown schematically in figure 5. 
The overall length of the model is 73.6 in. and the 
inlet height and width at the cowl tip are 8 and 
6 in. Initial flow compression is obtained with 7° 
compression ramps on the top and bottom walls 
over the first 15.6 in. of the inlet section. This 
setup is followed by 3.5° sidewall compression ramps, 
starting 21.6 in. from the inlet mouth. The sidewall 
compression extends 10.8 in. to the isolator section 
of length, 17.2 in. 

Model Configuration and Test Conditions 

The facility was configured to simulate the inlet 
flow conditions for a hydrogen-burning scramjet be- 
hind the bow shock of a 12° half- angle conical body 
immersed in a Mach 4 free-stream flow. Nominal val- 
ues were 1640°R for the total temperature, 92 psia 
for the total pressure, 11 psia for the dynamic pres- 
sure, and a Mach number of 3.5 for the gas at the 
model inlet. 


The model was mounted such that the inside top 
wall surface was continuous with the inside top wall 
surface of the square tunnel nozzle. This approach 
permitted simulation of the boundary layer growth 
along the upper inlet wall of an integrated engine 
module. At the isolator entrance the flow is com- 
pressed such that the core Mach number is reduced 
from 3.5 atjbe inlet to about 1.8. The flow from the 
isolator section enters the combustor section where 
fuel struts are located at the area step. As indi- 
cated by the perspective diagram in figure 4, these 
struts can provide a combination of tangential and 
transverse injection of gaseous hydrogen; however, 
only transverse injection was used in this test. Also 
not used in this test were an internal splitter and a 
fuel strut upstream of the area step in the isolator 
section shown in figure 4. During combustion, the 
heated flow expands through the 10° nozzle into the 
exhaust duct. To measure performance, the model 
was equipped with 16 calorimeter ports, 204 static 
pressure taps, and flowmeters to measure internal 
flow and fuel consumption. A standard wind tunnel 
balance was used to measure thrust. 

Pressure Gage Installation 

High-frequency fluctuating pressure measure- 
ments beneath high-enthalpy boundary layer flows 
present a challenging problem. Commercially avail- 
able high-frequency gages were employed that used 
piezoresist ive transducing elements. The gages were 
temperature compensated from 80° to 450°F, with 
an upper temperature limit of about 600° F. To re- 
duce sensitivity changes and probable destruction 
from thermal stress, it w T as necessary to either ac- 
tively cool the gages or recess them into the wall. For 
mechanical convenience, recess mounting was chosen 
in this test, with due consideration for the possibility 
of spurious noise and distorted frequency response. 

During a selected number of model performance 
tests, nine calorimeter port locations were made 
available for piezoresistive pressure gages distributed 
along the entire length of the model plus one gage 
in the lip of the tunnel nozzle. Six of these gage 
installations are shown in the distorted-scale side 
view of figure 6. The 0.159-in. gages were recessed 
from the interior wall surface by amounts varying 
from 0.78 to 1.97 in. to minimize damage from in- 
cident thermal flux to the diaphragms. The recess 
hole diameters were 0.170 in. Although gages were 
installed in the combustor and exhaust regions, lim- 
ited gage life prevented the acquisition of meaningful 
data. However, the gage at location 6, in the for- 
ward region of the combustor, did survive. Gage de- 
struction apparently was caused by hot gas spillage 
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flow around the outside of the model during un- 
starts, as opposed to excessive heat flux incident from 
the boundary layer. The surviving gages were lo- 
cated relative to the model inlet entrance at distances 
of —2.2, 8.5, 24.2, 35.0, 47.7, and 53.4 in. at recess 
depths of 0.78, 1.13, 1.97, 1.97, 1.85, and 1.75 in., 
respectively. All gages were designed for a pressure 
range of 0 to 100 psia. One other piezoresistive gage, 
with a range of 0 to 1000 psia, was installed down- 
stream from a quick-acting valve in a parallel branch 
of the fuel supply line. This gage provided a fuel 
pressure time history for correlating changes in com- 
bustion noise to incremental changes in the fuel flow 
rate. 

Additional tests were performed in a Mach 4 
blowdown tunnel (cold flow) to explore sensitivity 
changes during thermal transients and recession ef- 
fects on gage response. Thermal transient effects are 
addressed in appendix A. Gage recession effects on 
the measurement accuracy are discussed in appen- 
dix B, which also contains a correction procedure 
that accounts for effects of flow interaction with the 
open ports of the gage-cavity system. Finally, possi- 
ble flow noise from the interacting boundary layer 
and port-cavity system (self-noise) is discussed in 
appendix C. 

Data Acquisition, Reduction, and 
Analysis 

An important objective of this test was to cor- 
relate transient events in the pressure time history 
at various gage locations with changes in the overall 
engine operating parameters. The use of 0- to 100- 
psia pressure gages allowed fluctuating pressure data, 
superimposed on the mean static pressure, to be ac- 
quired at each gage location. In this manner the 
transient events could be visually analyzed and cor- 
related in context with the local mean static pressure 
at each gage location. 

The data were reduced in three stages. First, 
time domain analysis focused on pressure transients 
associated with engine operation transients such as 
onset of fuel flow, ignition, inlet unstart, and sudden 
heat release. Second, pressure spectra were obtained 
on pressure time history segments of special interest. 
Third, probability density functions were calculated 
on these same data segments. In this section, the 
acquisition, reduction, and analysis procedures will 
be discussed. 

Data Acquisition System 

Figure 7 is a diagram of the instrumentation used 
to obtain the fluctuating pressure measurements. 


Gage excitation and signal conditioning were sup- 
plied by a precision, low-noise, signal conditioning 
system below the tunnel test cell. The low-level 
signals from the pressure gages were transmitted 
through about 50 ft of shielded cable to the signal 
conditioner. Signals from the signal conditioner sys- 
tem were transmitted through about 75 ft of shielded 
cable to the digital signal acquisition system in the 
control room above the test cell. Fixed gains for each 
channel could be preset at the signal conditioner by 
changing the resistors in a feedback loop. 

The digital data acquisition system consisted 
of programmable bandpass filters, transient data 
recorders (TDK’s), and a computer system. The 
bandpass filter, with a roll off of 70 dB per octave, 
was employed as an antialiasing filter; the filter also 
had programmable gain that optimized the discrete 
signal level. The TDR’s discretized the output signal 
levels from the filters ( — 10 to 10 V) into 65 536 steps. 
The discrete samples were then temporarily stored 
in 0.5-megasample/channel buffers associated with 
the TDR’s. Depending on total gain selections, 
the pressure resolution ranged from 6.1 x 10"” 5 
to 94.3 x 10 -5 psia per step, for a nominal intrinsic 
gage sensitivity of 0.71 mV/psia. The acquisition 
procedure was controlled by the computer. Data 
sampling rate, filter cutoff frequencies, and signal 
amplification were inserted into the program as con- 
trol parameters. A more complete description of the 
data acquisition system and analysis software is given 
in reference 4. 

During a typical test run of approximately 28 sec, 
the pressure time history would generally include a 
number of transients separated by segments of statis- 
tically stationary pressure fluctuations. To capture 
all the data during a test run, the sampling rate was 
set to 15 kHz (data acquisition duration of 33 sec) 
and the bandpass filter was set in a low-pass mode 
with the cutoff frequency chosen to permit an analy- 
sis bandwidth of 5 kHz. This configuration permitted 
events such as inlet start, fuel ignition, and inlet un- 
start to be captured. After each test run, the TDR 
buffers were downloaded into permanent storage and 
a quick-look analysis of the data was performed. 

System Calibration 

The frequency response of the gage cavity sys- 
tem was expected to differ significantly from that 
of a flush-mounted gage. As a first step in evaluat- 
ing these effects, the frequency response of the gage 
cavity system was measured with a small-diameter 
(0.2-in. by 0.4-in.) waveguide that incorporated a 
reference gage and a high-frequency sound source 
described in reference 5. This frequency response 
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measurement, however, did not account for flow ef- 
fects. Flow effects on the recessed gage frequency 
response were examined analytically. This issue is 
addressed further in appendix B. 

The pressure gages were calibrated in situ with 
a precision constant pressure source at various times 
during the tests. A two-point calibration was used 
that included the endpoints of the expected pres- 
sure range at a specific gage location. This direct 
current (dc) calibration was used to adjust the fixed 
gains at the signal conditioner to accommodate the 
expected range of pressure during a given test run. 
The midpoint of these pressure ranges (i.e., the mean 
static pressures) was estimated from expected de- 
sign performance of the engine model. Consequently, 
there were occasional out-of-range conditions at var- 
ious gage locations that could not be accommodated 
by the a priori chosen dynamic range. This condition 
generally occurred more frequently for gages farther 
from the inlet, particularly in the combustor region. 

Data Reduction 

Of the 64 test runs during which fluctuating pres- 
sure data were acquired, the tunnel operating con- 
ditions were nominally the same; only the engine 
operating conditions were varied. Many of these 
runs were discarded after a visual inspection for suf- 
ficiently long segments of nominally statistically sta- 
tionary data from which pressure spectra could be 
calculated. From the surviving test runs, three were 
chosen as representative and will be discussed in some 
detail. 

Pressure spectra were calculated using standard 
time-series analysis procedures described by Hardin 
(ref. 6) and implemented with the same software 
as was used in reference 4 (see that reference for a 
discussion of the specific software routines). 

Data Analysis 

An immediate goal of the data analysis was to 
correlate events in the pressure time histories at dif- 
ferent gage locations with inlet unstarts. This cor- 
relation was done by visual inspection of the pres- 
sure time histories and by measurement of shock 
wave propagation speeds between gage locations. 
Also, measurements of pcak-to-peak fluctuating pres- 
sures yielded rough estimates of the rms fluctuating 
pressure differences between gage locations. 

Pressure spectral estimates were obtained during 
those portions of the pressure time histories that, on 
a visual inspection basis, approached statistical sta- 
tionary. Because the gages were mounted in re- 
cessed cylindrical cavities, directly measured power 


spectra were subject to distortion from cavity reso- 
nances, flow interaction with the port, and self-noise 
contamination generated by the grazing flow over the 
cavity port. Thus, some attention was given to these 
effects. 

Tijdeman and Bergh (ref. 7) describe a correc- 
tion procedure whereby the distortion effects alluded 
to above were correctly modeled for subsonic flow 
and low frequencies. A modified version of that 
model was investigated for the high-speed flow, high- 
frequenc y c ase; this modification was based on a 
suggestion by Heller and Widnall (ref. 8). Appli- 
cation of Heller and Widnall’s modification to the 
Tijdeman and Bergh model for the range of aerother- 
mal flow parameters and excitation frequencies of in- 
terest in this investigation resulted in corrections of 
not more than ±1.5 dB. Because this correction was 
well within the random measurement error range, it 
was not applied to these results. Essentially, the 
model calculates the transfer function between the 
fluctuating pressure incident on the cavity port and 
that at the face of the recessed gage. It takes into 
account the damping effect of viscous losses on the 
wall of the cavity and the pressure drop caused by 
the interaction of the steady flow and the fluctuating 
flow into the cavity port. (See appendix B for a more 
thorough discussion of the model.) 

The effects of spatial averaging on the finite sens- 
ing area of a transducer (in this case the port area) 
will reduce its effective sensitivity with increasing fre- 
quency. Blake’s review of sensor spatial averaging ef- 
fects (ref. 9) suggests that for port diameters used in 
this investigation, loss of signal from spatial averag- 
ing was less than 1 dB for the highest frequency of 
interest (5 kHz) in this investigation. 

Power spectral densities (PSD’s) were calculated 
for selected portions of time histories with accept- 
able statistical stationarity. Acceptable stationarity 
is largely an engineering judgment based on a com- 
parison of the means and variances of data blocks 
composed of a segment of pressure time history sub- 
ject to PSD analysis. In fact, no pressure time his- 
tories were ideally stationary. However, every effort 
was made to limit the analysis region such that the 
trade off between statistical accuracy and nonstation- 
arity effects was a balanced one. To improve statis- 
tical precision, the statistically stationary portions 
of the time histories were subdivided into contiguous 
data blocks of 4096 samples. Thus, for a sampling 
rate of 15.6 kHz, the data blocks are 0.262 sec long. 

These blocks were very nearly statistically inde- 
pendent as determined by the first zero crossing of 
the autocorrelation function. Subdivision of a data 
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record into Ng statistically independent data blocks, 
each of length tg, allows power spectral density esti- 
mates taken from each block to be treated as a chi- 
square random variable with two degrees of freedom. 
Thus, the larger Ng can be made, the smaller the 
statistical uncertainty. In other words, the accuracy 
of the power spectral estimate S p (lj) depends upon 
finding a sufficiently long stationary segment of data 
to provide a value of Ng that is as high as possible. 
Generally, Ng ranged from 7 to 31 blocks. Thus, the 
statistical accuracy of S p (<jj) was somewhat variable. 
With 7 blocks, the probability is 80 percent that 
S p (ij) lies between 0.66 and 1.72 of the true value; 
for 31 blocks the probability is 80 percent that the 
estimate lies between 0.81 and 1.30. The power spec- 
tral estimates were averaged over these Ng blocks as 
follows: 

1 Nb 

= (!) 

The trade off for reducing statistical uncertainty 
by increasing Ng is reduced frequency resolution 
bandwidth of the spectral estimate for a given total 
time length £ tota j. This trade off is quantified by the 
equation 

k = 2A / < total (2) 

where k represents the statistical degrees of freedom 
and A / is the resolution bandwidth associated with 
the data block length tg (A/ = 1 /tg). 

The discrete Fourier transform estimate for each 
data block was computed from 

A t N ~ l 

S p(u) ~ 2tt T) p (™) exp (~ iujn At) ( 3 ) 

n = 0 

where N — 4096. In this equation a “boxcar type” 
data window has been assumed that is unity over 
tg and zero elsewhere. Thus, the final spectral es- 
timate obtained by averaging the individual spectral 
estimates from Ng data blocks is given by 

S p (w) = W s |S p (u;)| 2 (4) 

where W s is a correction factor that accounts for the 
data window and experimental calibration factors. 
The rms pressures were computed by an integration 
of the area under S p (u;) from 1.9 Hz to 5 kHz. (The 
contribution of the frequency band centered at 0 Hz 
was excluded to avoid nonstationarities arising from 
slow variations of the mean static pressure.) 

Statistical Analysis 

In some instances, it is useful to compare mea- 
sured probability density functions with normal den- 


sity functions based on the computed mean and 
variance of the time history data as follows: 

nn b 

"-jfsjE'W (5) 

n— 1 

, AW. 

° 2 = nn ^ < 6 > 

Note that these statistical parameters are computed 
across all NNg values of the Ng contiguous data 
blocks of N data points each. Another measure of 
the rms pressure is the square root of the variance. 
However, this measure implicitly involves the mean 
as indicated in equation (6). Thermal transients have 
been shown to alter the mean values recorded with 
gages of the type used in this study; therefore, it was 
decided that an alternate form of this value should be 
calculated from the power spectral density with the 
contribution from the mean excluded. In all prob- 
ability density functions shown in the section “Dis- 
cussion of Results,” the tabulated overall sound pres- 
sure level (OASPL) was taken from the variance. In 
most cases, this value is close to the OASPL com- 
puted from the PSD (denoted (OASPL) psd ); how- 
ever, there are cases where the value is substantially 
different. This difference is an indication of a chang- 
ing underlying mean value (either “true” or ther- 
mally induced) or a lack of stationarity. When they 
arise, these situations will be noted. 

Skewness, which measures asymmetry in the 
density function, is calculated from 

^3-M3/^ 3 (7) 

where gs is the third moment of the pressure proba- 
bility density about the mean normalized by the third 
power of the standard deviation a, Kurtosis, which 
measures peakedness or flatness of the density func- 
tion, is calculated from the fourth moment normal- 
ized by the fourth power of the standard deviation; 
thus, 

a 4 = fi 4 Ja 4 (8) 

where g 4 is the fourth moment of the probability 
density about the mean. 

The calculated mean and variance were used to 
construct a Gaussian probability density function 
(PDF) against which the measured density function 
was compared. The measured density function is 
displayed as a normalized 25-bin histogram centered 
on the mean static pressure normalized by the rms 
pressure, as determined from the variance. 
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The section “Discussion of Results” contains an 
explanation of deviations in the PDF’s from Gauss- 
ian that is based upon several simulated time his- 
tories generated when a base or dominant pseudo- 
Gaussian process was perturbed by another whose 
mean is within the base distribution range but whose 
variance is much smaller. Such perturbations can be 
studied by observing changes in skewness and kur- 
tosis when the mean and variance of the perturbing 
process are altered. A high kurtosis value, relative 
to 3 for a perfect Gaussian, indicates that the instan- 
taneous pressure remains for a relatively long period 
near the mean value; a nonzero skewness value, rel- 
ative to zero for a perfect Gaussian, indicates that 
asymmetric excursions from the mean ^ccur. Fgr 
the experienced observer, these statistical parame- 
ters offer useful ways to monitor changes in the fluc- 
tuating boundary layer pressures that are not readily 
discernible in the power spectra. 

Estimates of Aerothermodynamic 

Parameters 

Aerothermodynamic parameters such as Mach 
number, unit Reynolds number, static pressure, and 
temperature were estimated at the various gage lo- 
cations from measurements of the tunnel operat- 
ing parameters at the model inlet and from static 
pressure measurements along the model wall. Esti- 
mates of aerothermodynamic parameters for the core 
flow at each measurement location were calculated 
on the basis of isentropic flow relationships. These 
estimates can err on the high side because of to- 
tal pressure losses through the developing oblique 
shock wave system. However, for these experiments, 
this error was not considered serious. For complete- 
ness, the isentropic flow relations for a thermally and 
calorically perfect gas, taken from reference 10 and 
appropriately manipulated for this work, are given 
below\ 

The Mach number M e at the boundary layer edge 
over a given gage location is given by 



where p e is equal to the measured static pressure p s 
at a wall port underneath the boundary layer and pt 
is the tunnel total pressure at the model inlet plane or 
tunnel nozzle exit (see fig. 6), The dynamic pressure 
q e is computed from 


Qe = 1 IPeMl (10) 


The static temperature at the boundary layer edge 
is given by 


T e = T t 



( 11 ) 


Unlike the static pressure, the static temperature in- 
creases through the boundary layer and will approach 
the adiabatic wall temperature for a nonconducting 
wall. In the present case the model mass was large 
and highly conducting, so thermal equilibrium was 
nev er reached and the wall temperature remained be- 
low the adiabatic value. The ratio of the adiabatic 
wall temperature to the static temperature in the 
core flow is given by 


1 -hi?/ 



( 12 ) 


For a nominal Mach number of 3.5 and a recovery 
factor Rj of about 0.9, an estimate of the adiabatic 
wall temperature for the flow through the present 
model is approximately 1460°R. 

Measurements of wall temperature indicated a 
mean of about 600° R. Thus, the ratio T w /T aw is 
about 0.4, which suggests a modestly cold wall con- 
dition. The gas density in the gage recess cavities is 
calculated from the gas law as follows: 


Pe 

pcav “ RT W 


(13) 


where the gas temperature in the cavity is taken to 
be T w but may be considerably higher because of gas 
exchange with the relatively hot boundary layer. 

Discussion of Results 

Results from three tests are presented here. These 
results are representative of the unsteady or fluctu- 
ating pressure loads in the GHSE for various operat- 
ing conditions. During the first test, run 6, batch 7 
(R6B7), a thermally induced unstart was captured. 
This test was interesting because it offered the op- 
portunities to observe an unstart-induced traveling 
shock wave over a sufficient distance and to de- 
termine the shock speed. However, only one gage 
remained within the measurement range after the un- 
start. Thus, a second test is presented that also in- 
volves an unstart; run 15, batch 6 (R15B6) has data 
records of sufficient length to provide high-quality 
PSD’s before and after unstart. This run was not as 
convenient for the determination of the shock speed 
because the surviving gages were too closely spaced. 
The third test, run 50, batch 10 (R50B10), was in- 
cluded to study the effects of step increases in the 
heat release rate without the complicating effects of 
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an unstart. This test is particularly interesting be- 
cause it provides evidence of combustion noise appar- 
ently propagating upstream in the subsonic portion 
of the boundary layer underneath the supersonic core 
flow in the isolator section. 

Table I lists information relevant to the test con- 
ditions and data analysis results. For each test run, 
the column entries were divided into four groups; the 
first vertical column of each group is the gage lo- 
cation number. The second vertical column of the 
first group is the model calorimeter port designation 
and nearest electronic scanning pressure (ESP) static 
pressure port. This information is useful to readers 
who want to correlate results from this report with 
those in the generic high-speed engine (GHSE) per- 
formance report of reference 3. The third vertical 
column presents gage location relative to the inlet 
leading edge, and the fourth vertical column gives 
the gage recess depth. The fifth, sixth, and seventh 
vertical columns relate to the specific segment of the 
measured pressure time history that was analyzed. 

The second horizontal group of entries contains 
the tabulated values of the thermodynamic param- 
eters and is generally self-explanatory with the aid 
of the symbols list. Vertical columns pt and 7} are 
stagnation values of pressure and temperature at the 
model inlet plane as estimated from tunnel oper- 
ating parameters. The next vertical column is the 
static pressure p e at the boundary layer edge, as de- 
termined from static pressure measurements p s at 
the wall static pressure ports. These values were ac- 
quired from averaging measurements acquired during 
the indicated time region as listed in the first table 
grouping. The next two vertical columns, labeled T e 
and are computed from the previous columns as 
discussed in the section “Data Acquisition, Reduc- 
tion, and Analysis.” The last vertical column, p c av , 
is computed from the gas law, with the temperature 
in the cavity taken to be 600° R as measured by a 
thermocouple installed in a calorimeter port. 

The third and fourth horizontal groups of en- 
tries relate to the correlation of unsteady pressure 
loads data with data for other fluctuating bound- 
ary layer loads in the literature. These entries 
will be referenced as further discussion develops. 
All dimensional quantities are expressed in British 
gravitational units. 

Run 6, Batch 7 (R6B7) 

Figure 8 shows curves for five axial static pressure 
distributions along the interior top wall of the model 
for test R6B7 just before and during an inlet unstart 
that is induced by thermal choking. The table in that 


figure lists the time coordinate for each pressure dis- 
tribution, the corresponding ESP cycle, and equiva- 
lence ratio <f) (i.e., the fuel-to-air ratio divided by the 
stoichiometric ratio). The column labeled cycle has 
the same information as the sixth vertical column of 
the first horizontal group in table I. This information 
allows the data in the present report to be correlated 
with the engine performance data of reference 3. 

The first three static pressure distributions, in- 
dicated by series 1, 2, and 3, correspond to a nor- 
mally operating engine (i.e., a net thrust production 
at an approximate equivalence ratio of 0.41). The 
static pressure rise is roughly linear with distance 
downstream from the model inlet to about 40 in. 
At that point, the pressure rise becomes somewhat 
steeper. This steeper region begins about midway 
in the isolator section. Maximum pressure occurs 
in the combustor region and decreases through the 
nozzle section. Note that mean pressure changes of 
2 to 3 psia (about 1.3 percent) occur in the combus- 
tor region over about 2-sec intervals; these changes 
are consistent with variations of about 2 percent in <fi. 

For curves 4 and 5, which correspond to equiv- 
alence ratios of 0.492 and 0.581, respectively, the 
pressure distribution is dramatically different. These 
curves reflect the unstarted state of the model. The 
onset of inlet unstart is associated with increased 
heat release rates as indicated by the higher equiv- 
alence ratios (see series 4 and 5 in the figure). The 
value of (j) for curve 4 is quite soft because the facility 
data acquisition system response time is not sufficient 
to update values of (f> in real time for this transient 
condition. For this particular test run, the equiv- 
alence ratio threshold for unstart is approximately 
0.58. In figures 9 to 14, the analysis is presented of 
the unsteady pressure loads data along the top wall 
of the model before and just after unstart. 

Figure 9 shows pressure time histories in 
pounds/square inch at gage locations 1 and 4 for 
test R6B7. Note that the use of absolute, as op- 
posed to differential, gages allows the high-frequency 
unsteady pressures to be viewed in context with the 
relatively slow mean pressure variations. The dif- 
ferences between the mean static pressures of fig- 
ures 8 and 9 (about 3 psia) at gage location 4 are 
likely due to temperature effects on the zero shift of 
the piezoresistive gages as discussed in appendix A. 
These effects notwithstanding, on the basis of argu- 
ments presented in appendix A, the static pressure 
variations in general were faithfully reproduced. 

The pressure fluctuations in figure 9 change 
throughout the test run, and significant events are 
indicated in the figure. Pressure fluctuations during 
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the initial 2 sec are indicative of tunnel and model 
inlet start transients and are not considered further. 
At 7 sec, the fuel is initiated, followed by ignition 
and thrust. The key feature, an inlet unstart, occurs 
at 13 sec, as shown by the sudden mean pressure in- 
crease at gage location 4 and is followed by a similar 
increase at gage location 1 exactly 0.436 sec later. 
As the disturbance emerges from the inlet, the signal 
conditioner for gage location 1 is driven out of range. 
These pressure increases indicate an upstream prop- 
agating shock that passes out the model inlet. The 
average upstream propagation speed of the shock be- 
tween the two gages, relative to the model reference 
frame, is about 7.1 ft /sec. As is evident at gage lo- 
cation 4, the fluctuating pressures increase dramati- 
cally after the inlet unstart. The fuel is turned off at 
17 sec and the facility is shut down at about 20 sec. 

A few points should be made concerning the mea- 
sured unstart phenomenon. First, the relatively slow 
speed at which the shock moves through the inlet in- 
dicates a potential for using control mechanisms to 
reverse the direction of the moving shock wave be- 
fore the unstart process is complete. Also, the conse- 
quence of an unstart is not just the increase in static 
pressure, but is also the exposure of an increased area 
of the engine to high fluctuating loads. 

If a researcher wishes to estimate power spectra 
and other statistical parameters of interest, it will 
not be possible to satisfy the statistical stationarity 
requirement for arbitrarily selected portions of the 
pressure time histories as shown in figure 9. There- 
fore, regions of the pressure time histories were cho- 
sen on the bases of interest and relatively good sta- 
tionarity for a sufficient length to provide meaningful 
averages of the PSD estimates over several contigu- 
ous data blocks. Two or three regions were. chosen 
for all three runs and are indicated by shaded regions 
in the pressure time history plots. Specifically, rel- 
atively good statistical stationarity is observed over 
intervals I and II from 3.7 to 5.5 sec and from 10 
to 12.3 sec, respectively. These time interval regions 
arc further analyzed. All aerothermodynamic entries 
in table I are averages taken over these regions. 

Figure 10 shows a magnified segment of the time 
history, 32 msec long, in region I at gage location 4. 
If the pressure fluctuations arc described approx- 
imately by a Gaussian process, then pressure ex- 
cursions remain within 3 standard deviations of the 
mean more than 99 percent of the time. There- 
fore, estimates of the rms pressures can be obtained 
by taking one-sixth of the peak-to-peak pressure 
fluctuations. Thus, a conservative estimate that 
is based on the indicated 8-msec time slice indi- 
cates a peak-to-peak value of at least 0.4 psia or 


147 dB (re 20 //Pa). A more accurate value obtained 
from the PSD over the 1.8-sec region I yields an 
(OASPL)psd of 153.1 dB. Thus, the particular peak- 
to-peak estimate illustrated in the figure is conser- 
vative by about 6 dB. Similarly, after unstart and 
before fuel shut off, peak-to-peak pressures approx- 
imate 7.0 psia (172 dB). The (OASPL)psp for the 
2.5-sec region III is 178.6 dB. 

Part (a) of figures 11 and 12 presents pressure 
spectra in two different formats at gage locations 1 
and 4 and for regions I and II of the pressure time 
histories. Pressure spectra for time region III (in- 
let unstarted) are not included because the gage at 
location f was out of range. The same data arc 
then shown in part (b) of the figures nondimension- 
alized for the flow velocity at the boundary layer 
edge t/ e , the rms value of the fluctuating pressure 
Prmsj and the boundary layer displacement thick- 
ness 8* = 0.23 in., as estimated at the tunnel noz- 
zle exit. The rms pressure is obtained by a sum- 
mation of the contributions from each frequency bin 
A/, except for the one centered at 0 Hz. In this 
way the contribution from the mean value is ex- 
cluded. The result of this calculation is listed as 
(OASPL)pgp> in the appropriate figures. All the 
power spectra are presented in this format to allow 
convenient comparisons with typical presentations of 
flat-plate boundary layer data. 

Figure 11(a) shows PSD’s expressed in decibels 
at gage locations 1 and 4 (sec fig. 9) for region I. 
The PSD analysis extends to nearly 8 kHz to ver- 
ify antialiasing filter roll off beyond 5 kHz. Beyond 
about 6 kHz, signal conditioning (discretization) 
noise floors, which are different for each instrumenta- 
tion channel, mask any further filtering action. These 
data above the filter cutoff are included to aid in 
judging the quality of the data below 5 kHz. For this 
case, the data below 5 kHz are well above the respec- 
tive channel noise floors and are apparently uncon- 
taminated by aliasing effects. Because of the varied 
ways these spectra may be affected by other contami- 
nating influences associated with recessing the gages, 
the reader is urged to consult the appendixes of this 
report. 

In region I, the inlet is started, but no combustion 
occurs. Except for the maxima near zero frequency, 
the general characteristic of the spectrum is that ex- 
pected for a turbulent boundary layer. The maxima 
near zero frequency may indicate upstream distur- 
bances in the tunnel or low-frequency model vibra- 
tion. Low-frequency vibration, however, may be the 
less likely contributor. Also, a relatively sharp, local- 
ized peak at 0.9 kHz is noted at reference gage loca- 
tion 1 on the tunnel nozzle lip. A broad peak is barely 
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perceptible near 2.1 kHz at gage location 4. The pos- 
sibility that either of these peaks is caused by gage 
mount vibration or organ pipe resonances of the gage 
recess cavity is discounted by arguments presented 
in the appendixes. For example, the first organ pipe 
resonance would be expected at 4.6 kHz if the gas 
temperature in the recess cavity were 600°R. Also, 
the pressure transfer correction model, discussed in 
appendix B, predicts a nearly complete damping of 
such resonances. Also, at 5 kHz, the level at gage 
location 4 has fallen only about 7 dB relative to the 
broad peak at 2 kHz. Clearly, much of the high- 
frequency energy above 5 kHz has not been captured. 

Comparisons of the two spectra show that the 
level at reference location 4 ranges higher than that 
at reference location 1 by 15 to 30 dB, whereas the 
overall level at gage location 4 is 18.7 dB higher than 
at gage location 1. From table I, the static pressure 
increases 5.24 times between the two locations, and 
Prms / Pe increases from 1.13 to 1.86 percent. These 
two factors account for the 18.7-dB increase in the 
rms pressure at gage location 4 relative to that at 
location 1. The static pressure increase is associated 
with the flow compression. The increase in Pnns/Pe 
may be associated with the interacting shock and 
boundary layers. Combustion is not involved in these 
spectra because there is none in this time frame. 

Figure 11(b) shows the same spectra in nondimen- 
sional form. Because each spectrum is normalized 
by its own rms pressure, the large differences alluded 
to above are largely suppressed. In this plotting for- 
mat, the differences between the two spectra are more 
clearly delineated. Also, the log-log scale empha- 
sizes spectral detail at the lower frequency range at 
the expense of detail at the higher frequencies. For 
example, the peak at gage location 1 near 0.9 kHz 
in figure 11(a) shows up near a reduced frequency 
of 3 x 10 -2 in figure 11(b). The spectrum at gage 
location 4, 35 in. downstream from the inlet plane, is 
somewhat more uniform with frequency than is the 
case for the spectrum at reference gage location 1, 
which is near the tunnel nozzle exit. In fact, the two 
spectra cross each other near a reduced frequency of 
7 x 10 -3 and diverge above that frequency, with the 
spectrum at gage location 4 nearly 10 times greater at 
a reduced frequency of 7 x 10~ 2 , near the filter cutoff 
frequency. This change in the boundary layer fluctu- 
ating pressure spectrum may be associated with the 
flow compression through the oblique shock in the 
inlet region. 

In the literature dealing with surface pressure 
fluctuations underneath turbulent flows, the free- 
stream dynamic pressure (e.g., refs. 11 and 12) is 
sometimes used to normalize the dimensional power 


spectra instead of the rms pressure as was done 
here. Other investigators use the dynamic pressure 
at the boundary layer edge as the normalizing fac- 
tor (e.g., sec ref. 11). The rms pressure was used 
here because it achieved the best possible collapse 
of the spectra from different gage locations and per- 
mitted similarities and differences to be more easily 
compared. 

Refer to table I to change the normalization pa- 
rameters. For example, if p rms is replaced by q c in 
figure 11(b), the PSD scale is multiplied by 2 x 10 -6 
and 23 x 10 -6 , respectively, for gage locations 1 
and 4. Also, the spectrum at gage location 1, re- 
gion I will follow the trend of the normalized spec- 
trum measured on a flat-plate equilibrium boundary 
layer underneath supersonic flow. References 1 1 to 13 
show that flat-plate spectra contain significant en- 
ergy to reduced frequencies near 1 for a wide range 
of flow and boundary layer parameters. Therefore, 
the present measurements miss at least 1 decade of 
the high-frequency spectrum. The fact that the spec- 
trum at gage location 4, region I is 14.7 times higher 
(and might even be higher with the missed portion of 
the spectrum included) suggests that the boundary 
layer turbulence grows more than would be the case 
for an equilibrium boundary layer on a flat plate for 
similar flow conditions. This assumption is further 
supported when rms pressure is compared as a per- 
centage of either the mean static pressure or dynamic 
pressure at the boundary layer edge. Table I shows a 
significant jump to 1.86 and 2.87 percent of the static 
pressure at gage location 4 (0.48 and 0.78 percent q e ). 

Figure 12(a) shows the pressure spectra at gage 
locations 1 and 4 in region II where combustion has 
reached steady state. The spectrum shape at refer- 
ence gage location 1 is practically identical to that 
for no combustion (region I) in figure 11(a). The 
(OASPL)psp> difference of only 0.7 dB between these 
regions suggests that very little combustion noise 
was measured at this gage location. At gage loca- 
tion 4, however, the spectrum, although quite similar 
in general shape and trend to that in figure 11(a), 
is up to 5 dB higher, starting well below 1.0 kHz 
and gradually coalescing at the high-frequency end 
near 5 kHz. Further, the (OASPL)p§p at gage loca- 
tion 4 in region II is 157.4 dB compared with 153.1 dB 
in region I (no combustion) at the same gage loca- 
tion and reflects a difference of 4.3 dB. This result 
suggests that broadband, low-frequency combustion 
noise may contribute to the fluctuating pressure at 
gage location 4. 

As was the case for region I, the gage location 4 
spectrum is about 30 dB greater than that at gage 
location 1 for the higher frequency range starting at 
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about 1 kHz. The spectral peak, near 0.9 kHz at 
gage location 1, is diminished somewhat; however, 
a peak at about 3.2 kHz becomes more prominent. 
These minor changes provide further evidence that 
the gages are measuring real pressure fluctuations 
rather than spurious output from gage mount vibra- 
tion. As with gage location 4 in figure 11(a), the 
pressure spectrum for gage location 4 in region II is 
devoid of the minor spectral peaks evident at gage 
location 1. 

Figure 12(b) shows the same data presented in 
figure 12(a), but in nondimensional form. The two 
spectra of this figure are also seen to collapse to 
a reduced frequency of about 7 x 10“ 3 , where a 
divergence begins and continues to the filter cutoff as 
was the case in figure 11(b). In this case, however, 
the divergence is more uniform and the broad peak at 
the reduced frequency of 8 x 10“ 2 is less noticeable. 
The pronounced effect of the filter roll off again 
suggests significant energy in these spectra beyond 
5 kHz. It is interesting to compare the spectrum at 
gage location 4 in this figure with the corresponding 
spectrum of figure 11(b) for the no-combustion case. 
There appears to be some increase in the level of the 
low-frequency spectrum at gage location 4, consistent 
with the increase in the decibel spectrum in region II 
(i.e., combustion). Also, the slightly less prominent 
peak near 10 -1 is probably a result of the increase in 
the gage location 4 spectrum between 10~ 2 and lO^ 1 . 
This finding is consistent with the comparisons on 
the decibel spectrum and suggests that broadband, 
low-frequency combustion noise affects the spectrum 
at the reference gage. 

Figures 13 and 14 show the probability density 
distributions for the pressure time histories at gage 
locations 1 and 4 and for time regions I and II. The 
continuous curves are Gaussian probability distribu- 
tions based on the means and variances as tabu- 
lated in the figures. The mean, variance, skewness, 
and kurtosis of the relevant regions have been cal- 
culated from the equations given in .sectio ns “Data 
Reduction” and “Data Analysis.” These values were 
calculated over all the data blocks used in the PSD 
analysis. The OASPL is calculated directly from the 
variance as listed. These values are also listed, along 
with those calculated from the PSD’s, in the last 
grouping of entries in table I. For the present test 
run, the quantity (OASPL — (OASPL)pgp) reaches 
a maximum of 3.9 dB, which occurs for time region II 
at gage location 1. For the most part, the differences 
between the two methods of calculating overall sound 
pressure levels are within 1 dB. The differences be- 
tween the two methods are likely caused by nonsta- 


tionarity in the pressure time histories from slowly 
varying mean static pressures. 

The measured probability distribution in fig- 
ures 13 and 14 is based on a 25-bin histogram cov- 
ering the pressure range. Each distribution has been 
normalized by its respective rms value (obtained from ; 
the variance), and the mean has been suppressed. 

In figure 13(a), the tabulated values for skew- 
ness and kurtosis of 0.15 and 2.96, respectively, sug- 
gest that the Gaussian density function describes 
the statistical properties of the pressure time history 
to a good approximation in region I. Visual inspec- 
tion of the data distribution for the Gaussian curve 
would suggest an inconsistency between the calcu- 
lated skewness and the apparent visual skewness. 
However, contributions to skewness are weighted by 
the cube of deviations from the mean. Thus, a careful » 
pairing of symmetrically occurring data points away j 

from the mean indicates relatively larger numbers of ; 

bin counts greater than 1 standard deviation from 
the mean p\ thus, p(t) — p/pvms = 1* 

However, bin counts beyond 3 standard devia- 
tions are of little significance because there are so few 
of them. Another item of interest in figure 13(a) is ; 
the tabulated OASPL of 135.2 dB. This value in no 
way depends on the Gaussian distribution assump- ; 
tion. It is computed directly from the variance of 1 
the measured time history values. The (OASPL)pgp 
obtained from the PSD was 134.4 dB. 

In figure 13(b), the PDF is shown for the same 
gage location but for region II of the time history 
where steady-state combustion occurs. Here the 
skewness is negligible, but the kurtosis has dropped 
to 2.70 (still 90 percent of that for a Gaussian curve). 

The visible data scatter about the Gaussian curve 
suggests a poorer fit than indicated by the skewness 
and kurtosis values. However, the apparent scatter 
is fairly sensitive to histogram bin width selection. 

For instance, if the bin width were greater, the ap- 
parent scatter distribution would change for adjacent 
values that are greatly different. The relation of vi- 
sual scatter to calculated values of kurtosis and skew- 
ness becomes increasingly difficult when the apparent 
data scatter becomes significant as in figure 12(b); 
however, if the data pairs away from the mean are 
compared, the symmetry appears to be controlled 
mainly by contributions between 1 and 3 standard 
deviations from the mean. The OASPL for time re- 
gion II is 139.0 dB as calculated from the variance 
(135.1 from PSD). This difference very likely reflects 
a lack of statistical stationarity. The PSD-based 
value is probably more accurate. 
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At gage location 4 (figs. 14(a) and 14(b)), the 
PDF’s are much more peaked (i.e., kurtosis greater 
than 3) and more skewed than at gage location 1. 
Therefore, the pressure excursions are more clustered 
near the mean value, as is evident in the figures. In 
regions I and II, the tails of the PDF’s show pres- 
sure excursions ranging from 3 standard deviations 
below to 10 standard deviations above the mean in 
region I. Even a few significant excursions from the 
mean pressure, particularly when an engine is oper- 
ating at nearly optimal conditions, may initiate the 
unstart condition. Thus, although the data for gage 
location 4 appeared mostly Gaussian, the excursions 
of approximately 10 standard deviations above the 
mean may provide insight into the stability of the air- 
flow through the engine (i.e., resistance to unstart). 

In region II, the skewness is about 3.7 times 
greater than in region I and the kurtosis is about 
2.7 times greater. These comparisons indicate that 
the pressure time histories become less Gaussian 
away from the inlet plane and the deviation from 
Gaussian becomes even more apparent when com- 
bustion occurs, in spite of the fact that the measure- 
ments are well upstream of the combustor and the 
core flow is supersonic. (Note also the listed overall 
levels of 153.2 and 157.8 dB calculated from the re- 
spective variances compared with 153.1 and 167.4 dB 
obtained from the corresponding pressure spectra 
previously described. The difference of nearly 10 dB 
in region II again suggests significant nonstationar- 
ity for the pressure time history in this region.) This 
situation is also consistent with the corresponding 
spectra of figures 11(a) and 12(a). This result sug- 
gests that combustion noise contributes significantly 
even in the inlet section of the model. 

Run 15, Batch 6 (R15B6) 

Run 15, batch 6 also involved an engine unstart. 
In this case, data were obtained at gage locations 4, 5, 
and 6. The gage at location 6, in the combustor re- 
gion (see fig. 6 and table I), survived long enough 
to acquire meaningful data during the combustion 
process. Figure 15 shows four static pressure dis- 
tributions along the model top wall before and after 
unstart in a similar format to figure 8. Also, the ESP 
system obtained a sequence of pressure distributions 
less than 1 sec apart that shows the progressive static 
pressure distribution changes during the unstart pro- 
cess. Series 1 and 2 are representative of static pres- 
sure distributions during normal model operation at 
equivalence ratios of 0.414 and 0.422, respectively. 
As the equivalence ratio approaches 0.5, the pressure 
changes as it did in figure 8, but with substantially 
different details. 


Figure 16 shows the pressure time histories at 
gage locations 4, 5, and 6 arid significant events in 
the test run. In this test, a third time history region 
is examined. This region occurs after the unstart 
event. Also, gage location 6 was out of range for 
pressure fluctuations corresponding to region I. 

Figures 17(a) and 17(b) show pressure spectra 
comparisons, expressed in decibels and nondimen- 
sional form, respectively, at gage locations 4 and 5 in 
region I before combustion was initiated. These spec- 
tra show similar trends, with the levels at gage loca- 
tion 5 slightly higher (fig. 17(a)) except near 2 kHz, 
where the gage location 4 spectrum peaks above that 
at gage location 5. Gage locations 4 and 5 were re- 
cessed to the same depth and the mean flow condi- 
tions were nearly the same over both gages. There- 
fore, any masking effects due to self-noise should be 
manifested to about the same degree in both spec- 
tra. Thus, the observed differences in spectral detail 
strongly indicate that self-noise did not significantly 
contribute to these spectra. Consistent with the 
trend established between gage locations 1 and 4, the 
(OASPL)psd increases from 153.4 to 156.8 dB. The 
nondimensional spectra again show similar trends, 
with gage location 4 levels exceeding those at gage 
location 5 to a reduced frequency of about 10~ 2 . 

Figure 18 shows a comparison of spectra at gage 
locations 4, 5, and 6 for region II during combustion. 
Combustion apparently causes the spectra at gage 
locations 4 and 5 to increase ((OASPL)pgjy y of 155.1 
and 165.2 dB, respectively); the greater increase 
occurs at gage location 5, with a consequent distinct 
separation between the spectra at gage locations 4 
and 5 across the entire frequency range. Also, except 
for the region below 2 kHz, the spectrum at gage 
location 5 is similar in trend to that at gage location 5 
before combustion. The core flow is supersonic in 
isolator section gage locations 4 and 5; therefore, 
combustion noise apparently propagates upstream 
through the subsonic part of the boundary layer 
to reach gage locations 4 and 5. More data are 
presented below to further support this contention. 
At gage location 6, in the combustor, the spectrum is 
characterized by five tone-like peaks to about 1.5 kHz 
and by a broad peak at about 2.4 kHz. Above 2.4 kHz 
the levels are generally 10 to 20 dB higher than at 
gage location 5 across the entire spectrum, with an 
(OASPL)psd of nearly 180 dB. 

Figure 18(b) shows the same data in nondimen- 
sional format. The spectra for gage locations 4 and 5 
collapse about the same as in test R6B7. However, 
the spectrum for gage location 6 deviates significantly 
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from the trends for gage locations 4 and 5; this re- 
sponse suggests that combustion noise in fact signifi- 
cantly contributes to this deviation as well as to that 
between gage locations 4 and 5. 

Figure 19 shows spectra at gage locations 4, 5, 
and 6 taken in region III of the pressure time his- 
tory after the inlet unstarted but before combus- 
tion ceased. After an unstart, the flow was subsonic 
throughout the model interior. The interesting fea- 
ture here is the high spectrum levels at gage loca- 
tion 4 ((OASPL)psd 177.5 dB and most distant 
from the combustor) relative to those at gage loca- 
tion 5 ((OASPL)psd of 167.5 dB). This phenomenon 
may be generated by acoustic coupling of the com- 
bustion process with longitudinal acoustic modes in 
the engine model. The standing waves resulting from 
such resonances may cause localized increases in the 
fluctuating pressure levels. Also, the spectrum for 
gage location 4 does not start to roll off until about 
5.8 kHz, which suggests a large spectral peak in this 
region that the low-pass filter did not attenuate suf- 
ficiently to avoid some aliasing back into the region 
to about 4.2 kHz. Thus, although the accuracy of 
the spectra for gage location 4 in this figure is cer- 
tainly compromised, the test results do suggest sig- 
nificant energy beyond 5 kHz; therefore, the over- 
all levels stated in this report are conservative. The 
same data are presented in nondimensional form in 
figure 19(b). Again, differences in these spectra are 
likely associated with combustion noise as previously 
discussed. 

Figures 20 through 22 show PDF’s for gage loca- 
tions 4, 5, and 6 for regions I, II, and III in the same 
format as for the previous test run. These PDF s ex- 
hibited the same general responses as did those for 
test R6B7 and will not be discussed in individual de- 
tail. It should be noted, however, that the highest 
rms sound pressure level of 180.3 dB was recorded at 
gage location 6 with low heat release (see fig. 22(a)). 
Also, a level of 178 dB was achieved at gage loca- 
tions 4 and 6 in the unstart condition (region III), As 
previously stated, these levels are almost certainly 
conservative in that significant contributions from 
frequencies above 5 kHz were filtered out at the data 
acquisition stage. Further, even though an acousti- 
cian w r ould view the fluid dynamic processes as quite 
nonlinear, the PDF s at gage location 6 in the com- 
bustor, both during normal engine operation with 
combustion and in the unstart condition, are well 
approximated by a Gaussian process (see fig. 22(b)). 
It should also be noted that the magnitude of the 
excursions from the mean was similar for all three 
gages and time regions. 


Run 50, Batch 10 (R50B10) 

The objective of run 50, batch 10 was to measure 
changes in fluctuating pressure levels from a sudden 
step change in the fuel flow and, thus, the heat re- 
lease rate. In this test an inlet unstart was avoided. 
Figure 23 shows the top wall static pressure distri- 
butions before, during, and after the step increase in 
heat release rate. Figure 24 shows a comparison of 
pressure time histories at gage locations 2 through 6, 
along with the pressure time history for a transducer 
installed in the fuel flow line. The time scale of the 
ESP was not completely synchronized to the pres- 
sure recording system that logged fluctuating pres- 
sures. Also, the gage at location 2 was mounted in 
an “infinite-line” configuration, similar to that pre- 
sented in reference 14. This configuration v*as chosen 
to provide additional thermal protection. 

Figures 25 through 27 show comparisons of the 
pressure spectra at gage locations 2 through 5 for re- 
gions I (no heat release), II (low heat release), and 
III (high heat release). These figures give evidence of 
combustion noise propagating upstream through the 
boundary layer. 4n particular, the pressure spectra 
shown in figure 27 show peaks at gage locations 4, 5, 
and 6 that gradually increase in frequency in the 
gages nearest the combustor region. Again, note 
that the change in gage recess depth is small com- 
pared with an acoustic wavelength and thus docs not 
sufficiently explain these shifts. 

Figures 28 through 32 show the PDF's at gage 
locations 2 through 6 for regions I, II, and III. Again, 
the same format was used as in previous PDF’s. The 
overall level for the high heat release region reaches 
181 dB (see fig. 31(c)) at gage location 5 as opposed 
to gage location 6 in the combustor. This somewhat 
surprising result is consistent with the (OASPL)psp 
computed from the PSD (171.5 dB and 167.7 dB at 
gage locations 5 and 6). It is interesting to note 
that the magnitudes of the excursions from the mean 
varied significantly versus location and heat release. 
This situation was especially true at gage location 2, 
where the scales of the plots had to be increased 
to include all excursions from the mean. Although 
some excursions were quite large, the engine never 
unstarted. This finding indicates that the engine 
was operating in a stable mode that was not easily 
disturbed. 

Table I summarizes all the key aerodynamic pa- 
rameters and rms fluctuating pressure data at the 
various gage locations. The entries in the table are 
useful for a comparison of the results of this inves- 
tigation with those from other investigations (e.g., 
refs. 11 through 13). The (OASPL) PS d and the 
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OASPL computed from the PDF variance generally 
agree well, as shown by the last three column en- 
tries for the last entry groupings. The agreement is 
particularly good for tests R6B7 and R15B7. For 
test R50B10, however, the differences are generally 
greater (to 12 dB), especially in high heat release re- 
gion III. These differences arc likely attributable to 
the relatively high degree of statistical nonstationar- 
ity in the pressure time history data in these regions, 
and that nonstationarity also affects the accuracy of 
the OASPL a s computed from the PDF. Thus, the 
(OASPL) psd computed from the PSD is likely more 
accurate. 

Table I also shows that the unit Reynolds num- 
ber increases from the inlet toward the combustor 
as the Mach number decreases. The ranges for the 
Mach and Reynolds numbers compare approximately 
with a portion of data in reference 12; in that ref- 
erence the Mach and unit Reynolds numbers ranged 
from 1.4 to 2.5 and from 7.2 x 10^ to 9.4 x 10^, respec- 
tively, to give rms pressures on a smooth tunnel wall 
from 0.2 to 0.4 percent q e compared with the present 
data for which rms pressure ranged from 0.13 per- 
cent q e at the model inlet to 7.5 percent q e in the 
isolator section. However, the rms pressure normal- 
ized by the mean static pressure ranged from 0.32 to 
nearly 17 percent, which is comparable to rms pres- 
sures of 5 to 24 percent of the mean static pressure re- 
ported in reference 13 for a flat plate in a hypersonic 
boundary layer flow ranging from Mach number 5.2 
to 10.4. Laganelli, Martellucci, and Shaw (ref. 11) 
show a comparison between a semiempirical predic- 
tion of rms pressures and measured data for attached 
turbulent boundary-layer flows for Mach numbers up 
to 10. Although much scatter is evident, both theory 
and experiment show a trend of decreasing p T ms/qe 
with increasing Mach number. In no case does mea- 
sured or predicted p rmii /q c exceed 1 percent. There- 
fore, the increasingly complex flow dynamics affects 
the rms pressure normalized by dynamic pressure 
which departs significantly from similar flat-plate 
data described in references 11 and 12. 

Figure 33 graphically summarizes all values of 
(OASPL) psd for the three tests and time regions 
in relation to gage location. At each gage location, 
the values of (OASPL) psd for the three time regions 
(I, II, III) are plotted from left to right. However, 
although a given time region designation is only 
a chronological ordering and does not consistently 
reference an engine operation state, these regions 
do include engine unstart conditions. The value of 
(OASPL) psd generally increases with gage location 
from the engine inlet. The data can be viewed as 
occupying a band of less than 30 dB at each gage 


location, ranging from 131 dB just upstream of the 
engine inlet to nearly 180 dB near the combustor. 

The last two figures (figs. 34 and 35) depict the 
skewness and kurtosis values versus axial distance 
(and gage location), again for all three tests and 
time regions. To keep the scaling acceptable, two 
kurtosis values were excluded (89.14 and 21.65 for 
test R50B10, gage location 2, regions II and III, 
respectively). Perhaps the most intriguing feature of 
these figures is that the data points were significantly 
non-Gaussian at gage locations 2 and 4, which were 
positioned at the start of both the inlet and the 
isolator. Thus, the non-Gaussian plot is possibly 
a reflection of changes in the shock wave structure 
where the wall angles change. Results of recent 
research (as yet unpublished) with a wedge on a flat 
plate indicate that this premise may indeed have 
some foundation. However, more data must be 
analyzed to determine the validity of this theory. 

Recommendations for Future 

Investigations 

The fluctuating wall pressure data in this report 
result from an exploratory investigation. Key is- 
sues were compromises in frequency response and dy- 
namic range associated with the recessed transduc- 
ers to achieve thermal protection. For frequencies to 
about 5 kHz, recession depths to about 2 in., and 
flow stagnation temperatures to about 1600°R, this 
approach appears possible. The gages were recessed 
in a massive heat-sink model with a ratio of about 0.4 
for the wall to adiabatic wall temperature. 

The high-speed flow over the gage recess ports 
may inhibit or dampen the organ pipe resonance of 
the recess cavities without generating excessive self- 
noise that could mask the desired fluctuating pres- 
sure measurements. However, measurements with 
both recessed and flush-mounted transducers should 
be performed to quantify any distorting effects of the 
open recess port on the spectrum, especially at the 
higher frequencies, or of small turbulence scales. 

Experiments also should be conducted to vali- 
date Tijdeman and Bergh’s transfer function (both 
amplitude and phase) with Heller and Widnall’s cor- 
rection for supersonic flows. That validation may 
make it possible to measure space-time correlations 
in high-temperature flows with recessed transducers. 

Conclusions 

This report has documented the results of an 
exploratory experiment to measure fluctuating wall 
pressures on an interior wall of a generic scramjet 
engine model. The tests described here were part 
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of a comprehensive series of performance tests con- 
ducted in the Langley Combustion-Heated Scramjet 
Test Facility, which can simulate Mach number 4 
flight conditions at altitudes ranging from 57000 
to 86000 ft. A key feature of the fluctuating pressure 
measurement approach was to recess the piezoresis 
tive pressure gages relative to the interior wall sur- 
face to protect them from the high-enthalpy flow. Al- 
though this procedure forced a careful interpretation 
of the data, the systematic error introduced by the 
recess apparently is not greater than about 1 dB. The 
measurements revealed several significant findings. 

1. The quality of the measurements suggests that 
recessed pressure transducers can be successfully 
employed to obtain meaningful fluctuating wall 
pressure data to at least 5 kHz in a hostile flow en- 
vironment with a tolerable transducer mortality. 

2. On the model wall near the inlet plane, root 
mean square (rms) pressures normalized by the 
local dynamic pressure at the boundary layer edge 
are comparable to those measured for attached, 
equilibrium boundary layers on flat plates for 
similar Mach and unit Reynolds number ranges. 

3. Pressure spectra generally are like those reported 
for flat plates near the model inlet, but are in- 
creasingly characterized by sharper peaks with in- 
creasing distance from the inlet plane. This pro- 
file is especially true near the combustor. These 
spectral peaks are apparently associated with the 
combustion process. 

4. With increasing distance from the inlet plane, 
including the isolator section, the rms pressures 
range to 7.8 percent of the dynamic pressure at 
the boundary layer edge for both normal engine 
model operation and the unstart condition. These 
fluctuating pressures are significantly greater than 
the upper limit of about 1 percent as reported 
in the literature on measurements for attached, 


turbulent boundary layers on flat plates for Mach 
numbers 0.4 to 10. Also, semiempirical the- 
oretical predictions for compressible, turbulent 
boundary layers over flat plates indicate decreas- 
ing Prms / Qe with increasing Mach number, with 
an upper limit approaching 1 percent at subsonic 
Mach numbers. These comparisons suggest that 
the fluid dynamics processes (such as interacting 
shock and turbulence and combustion noise) are 
more complex in producing scramjet fluctuating 
wall pressures than for flat-plate boundary layers. 

5. Comparison of pressure spectra at gages in the 
isolator section before and after increased heat 
release provides compelling evidence of combus- 
tion noise propagating upstream against the su- 
personic core flow. 

6. During the engine unstart process, the upstream 
shock wave propagates at about 7 ft/sec. There- 
fore, direct structural loading from pressure 
jumps associated with the unstart shock is likely 
benign relative to structural damage. 

7. On a decibel scale (re 20 /iPa), overall sound pres- 
sure levels range from 131 dB near the inlet plane 
to 180 dB just upstream of the combustor before 
an unstart condition. After an unstart, limited 
data suggest that a maximum level of 180 dB was 
not exceeded, although the spatial distribution of 
maximum levels may be different with higher lev- 
els extending into the isolator section, presumably 
because supersonic core flow was absent. 

All data presented in this report are subject 

to the systematic and random measurement errors 

discussed at length in the report and the appendixes. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
September 25, 1992 


16 


Appendix A 

Transducer Characteristics 

The transducer element for the high-frequency gages used to measure fluctuating pressure loads in this 
test was a piezoresistive transducer bonded to a metal diaphragm. The gage was packaged in a miniature, 
ruggedized case, capable of operating at up to 500°F. The gages were installed in the generic, high-speed 
engine (GHSE) with specially designed inserts to take advantage of existing calorimeter ports. A typical insert 
is illustrated in figure Al. The gage was mounted in a 0.159-in-diameter channel of the calorimeter insert 
as indicated in the figure. That insert was then flush mounted into the top wall of the GHSE model. This 
arrangement formed the recess gage mount used at all measurement locations in the GHSE model. The gages 
were recessed from 0.78 in. to 1.97 in. deep. Other relevant specifications for the gages are given below. 


Pressure range, psia 100 

Nominal sensitivity, mV/psia 0.75 

Combined nonlinearity, hysteresis, and repeatability at full-scale output, percent ... ±1 

Acceleration sensitivity: 

Perpendicular, percent full scale/g 0.0007 

Transverse, percent full scale/g 0.00014 

Compensated temperature range, °F 80 to 450 

Thermal zero shift, over compensated range, percent ±10 

Thermal sensitivity shift, over compensated range, percent ±7 

Diaphragm resonance, kHz 80 

Diaphragm diameter, in 0.147 


Because both the steady and fluctuating loads were to be measured, absolute gages were chosen instead of 
the differential type. Generally, the fluctuating pressures were expected to be at least an order of magnitude 
lower than the steady pressures. A two-point static pressure calibration, spanning the expected range of 
pressure excursions about the mean pressure at a given gage location, was applied to each gage. 

Also, thermal- and acceleration-induced sensitivity shifts were concerns. No vibration data were obtained 
during the tests. However, because of the massive GHSE model (37.3 slugs), any vibration-induced contamina- 
tion of the measured pressure fluctuations generally occurs at low frequencies. Based on the stated acceleration 
sensitivities, a worst case vibration of 100 Hz with an amplitude of 0.1 in. perpendicular to the gage diaphragm 
would generate a contaminating response equivalent to a 108-dB acoustic pressure. Generally, the low-frequency 
pressure spectrum levels recorded in the tests were at least 115 dB. Therefore, the contaminating influence of 
low-frequency vibration was not greater than about 1 dB. 

Comparison of mean static pressures measured by the piezoresistive pressure gages and by the facility ESP 
system during preliminary test runs indicated that the piezoresistive gages underwent substantial zero shifts, 
apparently due to transient thermal effects. Although the fluctuating pressure was measured for the general 
local mean static pressure, accuracy of the mean static pressure measurements was not considered crucial. 
There was concern, however, that measurement accuracy for the fluctuating pressure component would be 
adversely affected by thermal transients. 

To alleviate this concern as well as to compare recessed with flush-mounted gages, an aluminum, cold- 
flow scramjet model, 12 in. long, was modified to accept pairs of pressure gage inserts at three locations 
streamwise along a side wall as shown in figures A2 and A3. At each of the three locations a pair of gages 
could be aligned in streamwise or transverse direction to the flow. The inserts allowed the gages to be recess 
mounted to varying depths or be flush mounted. Exploratory tests were conducted on this model in the Mach 
4 Blowdown Facility to compare flush-mounted with recessed gage dynamic responses in a supersonic flow 
environment without thermal effects on a flush-mounted gage. Thermocouples were also installed in flush 
and recessed configurations similar to arrangements in the GHSE to provide information on gage exposure to 
thermal transients. 
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The Mach 4 Facility draws air from an outdoor tank farm that is some distance away. The tank farm 
air may vary daily by at least 20°F, depending upon meteorological conditions. Unfortunately, when the 
scramjet inlet was started in this facility, the boundary layer loads were not strong enough to generate coherent 
unsteady pressures that could be acquired by the present gages and analyzed. Therefore, the inlet was purposely 
unstarted to generate sufficiently high unsteady pressures for analysis (comparison of flush with recessed gages). 
This procedure provided some insight into the comparison of flush with recessed gage responses; however, the 
data were inconclusive because the aerodynamic environment was uncontrolled. Nonetheless, the gage responses 
to the mean static pressure were well matched to the static pressures recorded by the ESP system. 

Figure A4 shows the temperature changes at the flush and recessed thermocouples in degrees Rankine versus 
time in minutes. Tunnel startup is at 1.1 min and tunnel shutdown occurs off scale in this plot. Apparently, 
thermal equilibrium is achieved at about 3 min into the run. The initial temperature rise of about 10° R 
may be caused by higher temperature air being purged from the piping system before air from the tank farm 
is accessed. The initial temperature rise is followed by a steady decline for about 1.5 min, at which time 
the flush-mounted thermocouple reading rapidly drops, then recovers to about 10° R below the value for the 
recessed thermocouple. The recessed thermocouple apparently approaches its equilibrium value of about 497°R 
asymptotically. Generally, these thermal responses were expected except for the sudden transient disturbances 
at tunnel startup and for the flush-mounted thermocouple near the onset of the equilibrium temperature. After 
an initial precooling run, however, the model regulates to equilibrium temperature much more smoothly. 

Figure A5 shows a comparison of the pressure responses at a flush-mounted gage for a precooled with a 
nonprecooled model run. Note that the responses during tunnel startup transients are quite different, but 
in the equilibrium part of the run, the responses are about the same. This difference suggests that the gage 
temperature-compensating mechanism does perform after the thermal transients subside. To further validate 
the thermal transient effects on static pressure response, the tip of a gage was suddenly submerged in an 
ice-water bath under ambient pressure conditions. The results of two consecutive tests, starting with the 
gage at ambient temperature, are shown in figure A6. The two identical tests produced response transients 
that were quite different, but again, the equilibrium values approached the ambient pressure. These tests 
also suggest that the temperature-compensating mechanism for the piezoresistive gages does not provide for 
accurate pressure responses during a thermal transient. After the thermal transient has diminished and a 
new temperature equilibrium has been achieved, however, the pressure gages will again give accurate pressure 
responses. 

These tests suggest that both the zero shift A and the dynamic sensitivity B of the gage transfer 
characteristic may be affected during thermal transients as suggested in figure A7. Note that the parameters 
A and B are functions of temperature; for example, 

E = A(T) + B(T)p (Al) 

The linearized relationships between the fluctuating component of the gage voltage output E , the fluctuating 
temperature T, and the fluctuating pressure p were of interest here. Equation (Al) is thus linearized and solved 

for E to give _ _ 

E — (A^p -b Bpp)T + Bp (A2) 

where the parameters A and B arc implicit functions of time because of their possible temperature dependence. 
Clearly, the fluctuating pressure is contaminated by the thermal transient if the temperature derivatives Ap 
and Bp arc not zero. The temperature-compensating mechanism will ensure that the parameter B is held 
constant over the compensated temperature range after thermal transients have subsided. However, these 
qualitative tests suggest that the compensation does not hold Ap and Bp constant during thermal transients. 

It should be noted, however, that the current tests were conducted with the sensors mounted in a massive 
copper wall. Thus, thermal changes to the flow environment should slowly transfer to the gages. Therefore, 
although some measurements in this study may contain gradual shifts in their mean value, the fluctuating 
results should be unaffected. 
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Figure A2. Cold-flow scramjet model. 
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Figure A3. Pressure gage installation in cold-flow scramjet model. Linear dimensions are in inches. 



Figure A4. Thermal response of model with no tunnel precooling. 


21 







E 



23 



Appendix B 

Gage Recession Effects 

The transfer function between the fluctuating pressure at the inlet port of a recessed gage and that at the 
gage diaphragm is needed to correct the measured pressure spectra obtained with a recessed gage beneath a 
boundary layer flow, Tijdeman and Bergh (ref, 7) have developed a model for such a measurement configuration 
for the case when the core flow is subsonic and at low frequencies. The model is essentially a transmission 
line theory that accounts for the viscous dissipation in small tubes and the flow-port cavity interaction by an 
additional acoustic impedance. This modeling of the flow-port interaction is a greatly oversimplified description 
of an inherently nonlinear problem. However, Tijdeman and Bergh showed that the model was reasonably 
accurate for flow velocities to a Mach number of 0.8 and frequencies to 120 Hz. Richards (ref. 15) converted 
this theory to a more elegant matrix formulation. The intent in this appendix is to apply this transmission-line 
matrix formulation to recessed gage configurations and to extend the model to account for supersonic mean 
flow and/or high frequencies. 


Generalized Configuration 

Figure B1 depicts a typical recessed gage configuration for gage protection in harsh flow environments. For 
this configuration, the transfer function of interest is the ratio p g /ps- It. is assumed that the fluctuating pressure 
and volume velocity at the port entrance are known. According to reference 7, the boundary condition at the 
port inlet surface is given by 

[psl [1 CpsUe/s 

u s \ ~ [° 1 

where C is an empirical constant related to the effective port area, p s and U s are the fluctuating pressure and 
volume velocity just inside the port, and p s and U e are the static density at the port and mean flow velocity 
at the edge of the boundary layer over the port, respectively. The values of p\ and U\ are the fluctuating 
pressure and volume velocity just inside the port (location 1), and s is the cross-sectional area of the port. 
When the mean flow velocity is supersonic, Heller and Widnall (ref. 8) suggest that U e in equation (Bl) should 
be replaced with c, the adiabatic-phase velocity of sound. 

A matrix formulation can be used to describe the relation between fluctuating pressures and volume 
velocities at one end of a constant-area tube element and those at the other end as follows. 

/ p i+l \ ~ cosh Z Ci(w+ i) sinh (^j+i^i+l) 

\Ui + 1 / Z c ,(ti+1) sinh (^M+l^M+l) cosh (0M+ i*M+i) 

for i = 1,3, and 5. In this formulation, and are the complex propagation constant and length of 

element (z, i~\- 1). Z c is the tube characteristic acoustic impedance that depends on the element geometry. 

Thus, if the boundary conditions are known at the junction between the three constant-area elements and 
at the face of the gage, the desired transfer function {p 9 /p s ) can be determined. These boundary conditions 
are as follows: at the junction between the three elements, the pressure is constant and the volume velocity is 
continuous; for example, 

P2 = fa = P5 ( B3 ) 

U 2 -U 3 -U 5 = 0 (B4) 




The relation 



given in reference 15 is used to describe the conditions at gage location 6. In this equation, j — y/ 1, a; is the 
angular frequency, 7 is the ratio of specific heats, V is the volume of the cavity at the faee of the gage, cr is a 
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diaphragm deflection factor, and k is the polytropic constant (k = 1 for isothermal process, k = 7 for adiabatic 
process). Finally, because the cavity dimensions at the face of the gage are assumed to be very small relative 
to the shortest wavelength of interest, p g is assumed equivalent to pg. 

Recessed Gage Configuration 

Utilization of this transmission-line matrix formulation technique is convenient for modeling various high- 
frequency recessed pressure gage configurations in supersonic flow. The particular gage recess geometry in this 
study is depicted in figure B2. For this configuration, the derivation above can be simplified to a combination 
of equations (Bl) and (B2) and a boundary condition similar to that in equation (B5). 

Figure B3 depicts a comparison of measured data and model prediction (taken from ref. 8) for M = 0, 
M = 0.8, and low frequencies (to 120 Hz). As seen in the figure, the model performs quite well for those 
conditions. For the present test, however, it is important to understand the effects of sound and/or grazing 
flow over a recessed cavity at supersonic flow rates (to M = 3.5) and higher frequencies (to 5 kHz). 

A preliminary test was conducted to determine the effects of sound over a recessed cavity (no flow) for 
frequencies to 30 kHz. A recessed gage, with a recession length of 1.25 in. and a diameter of 0.16 in., was 
placed opposite a flush-mounted reference microphone in a small, rectangular waveguide (0.2-in. by 0.4-in.), in 
which only plane waves propagate over the frequency range of interest. The results from this test are shown in 
figure B4, where the solid curve represents the actual measured transfer function (ratio of acoustic pressure at 
the recessed gage to that measured by the flush-mounted high-quality condenser microphone) and the dashed 
curve represents the predicted transfer function. As can be seen in the figure, the model matches the measured 
levels at the three resonances. The frequencies where resonances occur are not quite as well matched but 
are remarkably close, given the sensitivity of the calculations. The model predicts the resonant frequencies 
much better than does the standard organ pipe calculation, which predicts that the resonances should occur 
at about 5, 15, and 25 kHz. Clearly, this transmission-line matrix formulation represents a useful model to 
account for various physical phenomena within the recessed gage cavity for higher frequencies (at least to 
30 kHz). 

Present Configuration 

Transfer function results are provided in figure B5 for a typical recessed gage configuration in the present 
test. The aeroacoustic parameters in the transfer function calculations are included in the inset table. The 
maximum frequency for the present test was 5 kHz, so the transfer functions are displayed only for that 
frequency range. In this figure, the solid curve represents a transfer function expected with no mean flow. The 
dashed curve depicts the transfer function computed for the supersonic flow condition. 

Note that the no-flow curve contains a resonance of 27 dB at approximately 2 kHz, whereas the supersonic 
flow curve indicates a peak resonance of less than 1 dB near 1.8 kHz. Because the transfer function for the 
supersonic flow case is so flat, this correction (recall that this transfer function typically should be used to 
correct the measured response to attain the true response) was not implemented in view of other greater 
experimental uncertainties. This situation held true for all port lengths in the experiment. 

Some measured spectra indicate peaks that fall near these predicted resonances. These measured peaks are 
often more than 1 dB above the surrounding levels, as predicted by the transmission line model (with all of 
the aeroacoustic parameters included), but are much less than 27 dB above the surrounding levels. However, 
no consistent trend appears for the resonant peaks in the measured spectra. Thus, this model may properly 
describe the effects of the recessed cavity, even for these harsh conditions. Theoretical and experimental studies 
are being conducted to improve the model for supersonic flow rates. 
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Figure B2. Recessed gage configuration. 
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Figure B3. 


Frequency response change of recessed gage from grazing flow (ref. 7). 



Figure B4. Measured versus predicted transfer function magnitudes and resonant frequencies 



Appendix C 


Recessed Gage Self-Noise Effects 

In this appendix, the self-noise or aeroacoustic disturbance generated by the interaction of the boundary 
layer with the gas in the recess port is addressed. Such interactions, if present, distort and contaminate the 
aeroacoustic loads in the absence of the recess port. The frequencies at which such effects may occur are 
estimated from results for subsonic flows over relatively shallow cavities. Reference 16 gives a survey of the 
work in this area. 

According to reference 16, flow-cavity interactions are classified into three categories: fluid dynamic, fluid 
resonant, and fluid elastic. Of these, the fluid-resonant interaction should be the dominant interaction mode 
for the recessed gage configurations in this investigation. (Recess diameter and length in this discussion are 
referred to as cavity length and depth, respectively, in the context of refs. 16 and 17. Also, the cavities of these 
references were generally rectangular.) Formulas quoted from reference 17 have been changed for consistency 
with the notation of this discussion, l is the recess (or tube) length, and d is the recess (or tube) diameter. For 
the recessed gages along the isolator section of the GHSE, the djl is about 0.08. 

Figure Cl shows a collection of data replicated from reference 17 for the first three crosswise modes (i.e., 
modes associated with the diameter d) for cavity djl ratios ranging from 16 to 0.25. The data are plotted as 
Strouhal number fd/U e versus Mach number M for three different mode numbers n . Although most of the 
data are for subsonic flows, a small amount extends to a Mach number of 1.1 for modes 1 and 2. The data arc 
fairly well described by an empirical equation attributed to Rossiter and given in reference 17 as 


fd n — a 

Ue = (W + M 


(Cl) 


The empirical constants K v and a are taken to be 0.57 and 0.25, respectively, for the first three modes. Note 
that the parameter djl does not explicitly appear in this equation. Reference 17 indicates that a decreases as 
djl decreases. The value K v is the ratio of vortex convection velocity across the recess port to the free-stream 
velocity. If this equation is used to calculate a Strouhal number for the flow conditions over a typical recessed 
gage of interest in the isolator section (i.e., M a* 1.6 and djl « 0.08), then for a first-mode excitation, the result 
is 0.22. This value appears to agree reasonably well with the extrapolation of the bottom curve in figure Cl to 
a Mach number of 1.6. This value for the Strouhal number gives a crosswise modal frequency of about 43 kHz, 
which will be higher if a lower value of a is used. This frequency is well out of our range of interest and is not 
a problem in this investigation. 

Flow excitation of depthwise modes in shallow cavities underneath subsonic flows is also discussed in 
reference 17; however, extrapolation to the configuration of interest in this investigation is far less reliable than 
is the case for the crosswise modes discussed above. In the context of the present investigation, these depthwise 
modes can be identified with the longitudinal acoustic or organ pipe modes that control the resonant response 
of the gage-cavity system to incident acoustic pressure disturbances discussed in appendix B. In addition to 
purely acoustic pressures, a possibility clearly exists for the unsteady part of the flow to excite these modes and 
to result in additional distortion of the desired fluctuating pressure measurements. However, as suggested in 
appendix B, localized mean flow interaction with the ported tube dampens the resonant response to negligible 
levels. Therefore, flow-port interaction noise is essentially uniformly distributed across the entire frequency 
spectrum; consequently, it just raises the effective noise floor of the measurement system. This flow-port 
interaction or self-noise is assumed present at about equal levels at all gage locations and is treated as the key 
contributor to the total noise floor of the measurement system. 
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Tabic I. Data for Recessed Gage Transfer Function Calculation 


(a) Run 6, batch 7 (R6B7) 


Model 


Port/ESP 

Location, 

in. 

TN/NOZ18 a 

-2.2 

TN/NOZ18 a 

-2.2 

306/51 

35.0 

306/51 

35.0 

vu 

Tu 

psia 

°R 

94.14 

1532.4 

92.38 

1580.2 

94.14 

1532.4 

92.38 

1580.2 


Analysis segment 


c ej 

ft /see 

1046.3 
1058.5 

1325.4 
1360.8 



U e , 

ft /sec 

3596.3 

3657.8 

3102.1 

3118.0 


2.65 x 10 6 
2.47 
5.22 
5.08 


Recession, 


(OASPL)psd> 

dB 


Region 

I 

II 

I 

II 


Pmis i 
psia 

0.0152 

.0165 

.1311 

.2150 

OASPL, 

dB 


ESP, 

Interval, 

cycles 

sec 

6-8 

3.7-5. 5 

15-18 

10.0-12.3 

6-8 

3.7-5. 5 

15-18 

10.0-12.3 


Pe? 

slugs/ft 3 

248.6 x 10“ 6 
232.1 
810.9 
815.4 


Pcavj 

slugs/ft 3 

188.8 x 10 
180.4 
988.3 
1047.6 


Prms/Pe j 

Prms/ Qei 

percent 

percent 

1.13 

0.14 

1.28 

.15 

1.86 

.48 

2.87 

.78 


OASPL - (OASPL)psd, 
dB 


NOZ18 refers to port 18 in tunnel wall. 




Table I. Continued 


(b) Run 15, batch 6 (R15B6) 


Model 

Port/ESP 

306/51 

306/51 

306/51 

308/81 

308/81 

308/81 

309/101 

309/101 


Analysis segment 


Location, 


Recession, 


Region 

I 

II 

III 
I 

II 

III 

II 

III 


ESP, 

Interval, 

cycles 

sec 

" 5-8 

3. 1-5.2 

14-20 

9.4-13.4 

24-36 

16.3-24.4 

5 8 

3. 1-5.2 

14-20 

9.4-13.4 

24-36 

16.3-24.4 

14-20 

9.4-13.4 

24-36 

16.3-24.4 


ft /sec 

1291.3 

1308.5 

1558.9 

1313.9 

1508.3 
1568.0 

1491.5 

1529.4 


1468.9 

1502.4 

1540.5 
1468.9 

1502.4 

1540.5 

1502.4 

1540.5 

U e , 

ft /sec 

3050.8 

3079.8 

2520.6 
3002.0 

2582.8 

2492.4 

2631.4 

2609.4 


5.29 x 10 6 
5.18 
7.27 
5.54 
7.26 
7.34 
7.09 
7.02 


(OASPL)psd, 

dB 


694.0 
712.7 

1011.6 

718.6 

947.0 
1023.4 

925.9 

973.7 


OASPL, 

dB 


Pe-> 

slugs/ft 3 

805.0 x 10 6 
795.9 

1737.9 

878.1 
1619.7 

1788.9 
1531.0 
1579.6 


Prms? 

psia 

Pnns/ Pc 3 

percent 

0.1357 

2.04 

.1650 

2.44 

2.1752 

10.38 

.2007 

2.67 

.5278 

2.89 

.6879 

3.15 

2.8674 

16.97 

2.3850 

13.01 


Pcavj 

slugs/ft 3 

931.1 x 10 
945.5 
2930.1 
1051.7 
2556.4 

3051.3 
2362.6 

2563.4 

P rms/ 9 ei 
percent 

0.52 

.63 

5.67 

.73 

1.41 

1.78 

7.79 
6.39 


OASPL - (OASPL) psd , 
dB 


Table I. Continued 


(c) Run 50, batch 10 (R50B10) 



Model 


Analysis segment 


Port/ESP 

301/7 
301/7 
301/7 
304/26 
304/26 
304/26 
306/51 
306/51 
306/51 
308/81 
308/81 
308/81 
309/101 



Location, 

in. 

875 
8.5 
8.5 
24.2 
24.2 
24.2 
35.0 
35.0 
35.0 
47.7 
47.7 
47.7 
53.4 


1475.7 
1575.0 

1581.7 

1475.7 
1575.0 

1581.7 

1475.7 
1575.0 

1581.7 

1475.7 
1575.0 

1581.7 
1581.7 


Recession, 

in. 

(5T 

(a) 

(a) 
1.97 
1.97 
1.97 
1.97 
1.97 
1.97 
1.85 
1.85 
1.85 
1.75 



Region 

T 

II 

III 

I 

II 

III 

I 

II 

III 

I 

II 

III 

III 


581.0 

623.1 
625.7 
694.0 

739.5 

742.5 

703.6 

770.9 

933.7 
866.6 

934.9 

937.2 
1050.9 


ESP, 

cycles 

Interval, 

sec 

4-6 

3. 1-5.0 

17-19 

14.2-16.0 

22 24 

18.3-20.2 

4-6 

3. 1-5.0 

17-19 

14.2-16.0 

22-24 

18.3-20.2 

4 6 

3. 1-5.0 

17 19 

14.2-16.0 

22 24 

18.3-20.2 

4 6 

3.1 5.0 

17-19 

14.2-16.0 

22-24 

18.3-20.2 

22-24 

18.3-20.2 




P cav, 
slugs/ ft' 5 


522.8 x 10" b 

506.3 x 10 

475.4 

493.7 

476.1 

496.5 

815.0 

942.7 

729.6 

899.3 

730.5 

904.0 

843.6 

989.3 

809.5 

1040.1 

1295.1 

2015.4 

1420.5 

2051.7 

1311.3 

2043.4 

1307.3 

2042.0 

1740.6 

3048.5 


Infinite line. 





Table I. Concluded 


(c) Concluded 



c e? 

E/e, 

Qe , 

Prms> 

Prms/ Pei 

Pxms/Qei 

Gage 

ft/sec 

ft /sec 

psia 

psia 

percent 

percent 

2 

1181.5 

3278.2 

19.5 

0.0115 

0.32 

0.06 

2 

1223.5 

3381.5 

18.9 

.0105 

.30 

.06 

2 

1226.0 

3388.7 

19.0 

.0109 

.31 

.06 

3 

1291.2 

3064.4 

26.6 

.0665 

.99 

.25 

3 

1332.9 

3167.9 

25.4 

.0781 

1.21 

.31 

3 

1335.6 

3174.9 

25.6 

.1341 

2.07 

.52 

4 

1300.1 

3045.4 

27.2 

.1708 

2.41 

.63 

4 

1360.9 

3107.9 

27.1 

.3369 

4.53 

1.24 

4 

1497.7 

2789.9 

35.0 

.9279 

6.44 

2.65 

5 

1442.9 

2704.8 

36.1 

.2007 

1.37 

.56 

5 

1498.7 

2772.8 

35.0 

.1388 

.95 

.40 

5 

1500.5 

2782.4 

35.1 

1.0902 

7.47 

3.10 

6 

1588.9 

2525.1 

38.5 

.6273 

2.88 

1.63 





N Iie > 

(OASPL)psd, 

OASPL, 

OASPL - (OASPL)psd, 

Gage 

M e 

ft" 1 

dB 

dB 

dB 

2 

2.8 

4.21 x 10 6 

132.0 

133.0 

1.0 

2 

2.8 

3.75 

131.2 

134.1 

1.9 

2 

2.8 

3.75 

131.5 

134.2 

2.7 

3 

2.4 

5.38 

147.2 

151.9 

4.7 

3 

2.4 

4.76 

148.6 

153.1 

4.5 

3 

2.4 

4.76 

153.3 

156.4 

3.1 

4 

2.3 

5.48 

155.4 

157.4 

2.0 

4 

2.3 

5.03 

161.3 

163.3 

2.0 

4 

1.9 

6.33 

170.1 

177.4 

7.3 

5 

1.9 

7.08 

156.8 

158.2 

1.4 

5 

1.8 

6.36 

153.6 

160.2 

6.6 

5 

1.8 

6.35 

171.5 

181.0 

9.5 

6 

1.6 

7.11 

166.7 

178.7 

12.0 
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Figure 3. GHSE model installed in Combustion- Heated Scramjet Test Facility. 
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Figure 4. Conceptual view of GHSE model. 
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Figure 7. Digital Data Acquisition System. 
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Series 


Time, sec 


Cycle 


Static 

pressure, 

psia 


1 

10.2 

15 

0.415 

2 

10.9 

16 

.407 

3 

12.2 

18 

.407 

4 

13.0 

19 

.492 

5 

14.3 

21 

.581 


a Acquired during transition to unstart. 



Figure 8. Axial static pressure distribution versus time for test R6B7. 


Gage location 1 




Figure 9. Pressure time histories during unstart for test R6B7. 
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Gage (OASPL)psd 



Frequency, kHz 


(a) Pressure spectra. 



(b) Nondimensional spectra. 

Figure 11. No-combustion case for test R6B7, Region I; 3.7 to 5.5 sec. 
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(a) Region I with no combustion. 
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(b) Region II with combustion. 

Figure 13. Probability density distribution of pressure time histories at gage location 1 for test R6B7. 
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(a) Region I with no combustion. 
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(b) Region II with combustion. 

Figure 14. Probability density distribution of pressure time histories at gage location 4 for test R6B7. 
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Figure 16. Pressure time histories during an unstart for test R15B6. 
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(a) Pressure spectra. 
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(b) Nondimensional spectra. 

Figure 19. Inlet unstarted case for test R15B6, Region III; 16.3 to 24.4 sec. 
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(a) Region I with no 

combustion. 

Gaussian 

Mean 7.00 psia 

Measured 

Variance .030 psia 1 


Skewness 0.58 


Kurtosis 6.85 

0 

OASPL 155.5 dB 


[p(t)-p]/ Pr , 


(b) Region II with combustion. 

Figure 20. Probability density distribution of pressure time histories at gage location 4 for test R15B6. 
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(a) Region I with no combustion. 
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(b) Region II with combustion. 

Figure 21. Probability density distribution of pressure time histories at gage location 5 for test R15B6. 
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(c) Region III with inlet unstarted. 
Figure 21. Concluded. 
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(a) Region II with combustion. 



Gaussian 
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(b) Region III with inlet unstarted. 

Figure 22. Probability density distribution of pressure time histories at gage location 6 for test R15B6. 
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(b) Gage locations 4, 5, and 6. 
Figure 24. Concluded. 
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(a) Region I with no heat release. 
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(b) Region II with low heat release. 

Figure 28. Probability density distribution of pressure time histories at gage location 2 for test R50B10. 
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(a) Region I with no heat release. 
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(b) Region II with low heat release. 

Figure 29. Probability density distribution of pressure time histories at gage location 3 for test R50B10. 
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(c) Region III with high heat release. 
Figure 29. Concluded. 
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(a) Region I with no heat release. 



(b) Region II with low heat release. 

Figure 30. Probability density distribution of pressure time histories at gage location 4 for test R50B10. 
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(c) Region III with high heat release. 
Figure 30. Concluded. 
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(a) Region I with no heat release. 



(b) Region II with low heat release. 

Figure 31. Probability density distribution of pressure time histories at gage location 5 for test R50B10. 
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Figure 32. Probability density distribution of pressure time histories at gage location 6 for test R50B10, 
Region III; High heat release. 
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Figure 34. Skewness versus axial distance for all tests. 
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Figure 35. Kurtosis versus axial distance for all tests. 
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